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SUMMARY 

In  a flexible  aircraft  flight  control,  load  control  and  structural 
mode  control  interfere  with  each  other.  Therefore,  an  integral  de- 
sign of  controller(s)  is  necessary.  This  paper  describes  how  an  in- 
tegral aircraft  model  covering  the  requirements  of  all  three  disci- 
plines can  be  derived  and  how  an  integral  controller  can  be 
designed  by  multiobjective  parameter  optimization.  General  de- 
sign criteria  for  mode  control  are  proposed. 

1.  INTRODUCTION 

Flight  control  and  aeroelastics  are  two  disciplines  for  aircraft  (a/c) 
design  which  have  worked  more  or  less  independently  from  each 
other  in  the  past.  Flight  control  deals  with  the  nonlinear  rigid- 
body  motion  of  the  a/c  and  aeroelastics  deals  with  linear  vibra- 
tions of  the  a/c  structure. 

This  was  possible  as  long  as  the  eigenfrequencies  of  rigid-body 
motion  and  structural  vibrations  were  clearly  separated  from  each 
other  and  slow  movements  of  control  surfaces  were  sufficient.  But 
with  modem  large  and  flexible  a/c  there  are  interferences  of  both 
types  of  motion. 

There  are  three  reasons  for  this. 

First  the  eigenfrequencies  of  rigid-body  motion  and  structural 
vibrations  come  close  together  causing  stronger  cross  coupling. 
Second  the  measured  input  signals  of  the  electronic  flight  control 
system  (EFCS)  might  contain  signal  parts  of  structural  vibrations. 
Feeding  back  these  signals  could  therefore  lead  to  instability  of 
flutter. 

Third  there  may  be  the  phenomenon  of  aircraft  pilot  coupling 
(APC)  where  vibrations  of  the  cockpit  floor  structure  may  lead  to 
movements  of  the  side  stick  via  pilot  seat  and  pilot  body /pilot  arm, 
thereby  causing  control  surface  movements  with  eigenfrequen- 
cies of  the  rigid-body  and  the  a/c  structure. 

With  modem  a/c  it  is  therefore  no  longer  possible  to  neglect  the 
coupling  between  rigid-body  motions  and  deflections  of  the 
structure.  Interdisciplinary  cooperation  of  flight  control,  aero- 
elastics and  of  the  discipline  loads  prediction  is  therefore  required. 
An  integral  a/c  model  and  an  integral  controller  covering  the 
requirements  of  all  these  disciplines  are  necessary  for  a successful 
a/c  design. 

In  the  following  some  details  of  integral  model  derivation  and 
some  remarks  on  the  design  of  integral  controllers  are  given.  Em- 
phasis is  put  on  the  structural  vibrations  part  of  this  controller. 

2.  THE  INTEGRAL  MODEL 

The  integral  mathematical  model  of  the  a/c  must  include  rigid- 
body  movements  and  structural  vibrations.  This  causes  problems. 
In  the  past  rigid-body  movements  studied  by  flight  mechanics 
covered  large  movements  such  as  large  angles  of  attack,  and  this 
required  nonlinear  aerodynamics  and  nonlinear  Euler-Newton 
equations.  Elasticity  or  structure  deformation  was  only 
introduced  via  ’’elastified”  aerodynamic  coefficients  taking  into 
consideration  steady  state  deflections  of  the  a/c  structure. 

On  the  contrary,  structural  vibrations  studied  by  aeroelastics  were 
described  by  linear  equations  paying  attention  to  the  first  harmon- 
ics of  unsteady  aerodynamics  with  small  amplitudes  of  angle  of 
attack  and  neglecting  completely  aerodynamic  forces  in  fore  and 
aft  direction.  Therefore,  rigid-body  motion  is  included  - if  at  all 
- in  these  linear  equations  only  in  a rather  rough  manner. 

An  integral  model  covering  the  requirements  of  both  disciplines 
may  therefore  be  rather  complicated.  It  should  superimpose  small 
deflections  of  the  grid  point  masses  of  a huge  FE-model  with 
large  three-dimensional  movements  in  space.  It  should  take  into 
account  aerodynamic  loads  due  to  large  angles  of  attack  and  large 
amplitudes  of  the  phygoid  for  100000  degrees  of  freedom.  In 
addition,  Navier-Stokes  CFD  codes  and  nonlinear  Euler-New- 
ton inertias  should  be  implemented.  But  for  the  time  being, 


another  approach  is  necessary. 

One  of  the  most  obvious  models  is  therefore  a linear  model  valid 
for  a working  point  under  steady  state  flight  conditions.  This  al- 
lows the  study  of  movements  around  this  point,  and  flights  within 
the  whole  envelope  are  possible  by  interpolation  between  differ- 
ent working  points. 

The  model  presented  here  is  based  on  the  traditional  aeroelastics 
model  in  frequency  domain.  But  it  has  been  improved  so  that  its 
rigid-body  modes  should  approach  the  modes  of  the  linearized 
flight  mechanics  model.  The  improved  model  includes  the  fol- 
lowing: 

-1  stiffness  matrix  of  structure 
-2  inertia  matrix  including  all  grid  points 
—3  aerodynamic  panels , also  for  the  fuselage 
— 4 steady  state  air  loads  at  the  working  point 
-5  steady  state  elastic  deformations  of  the  structure 
-6  corrections  of  the  stiffness  matrix  due  to  steady  state 
deformations 

-7  corrections  of  the  inertia  matrix  due  to  steady  state 
deformations 

-8  eigenvalue  analysis  of  the  structure 
( without  aerodynamics) 

-9  mode  reduction  by  truncation 
-10  addition  of  control  modes  for  compensation  of 
truncation  effects 

-1 1 transformation  of  rig  id-body  modes  from  ma in  inertia 
axes  to  geodetic  aircraft  reference  axes 
-12  traditional  unsteady  air  loads 
-13  addition  of  missing  elements  to  traditional  unsteady 
air  loads 

-14  inclusion  of  gust  loads 

-15  smoothing  of  rigid-body  air  loads  at  zero  frequency 

-16  engine  loads 

-17  weight  force 

-18  structural  damping 

-19  linearization  and  transformation  into  time  domain 
-20  eigenvalue  analysis  at  working  point 
-21  introduction  of  nonlinear  rigid-body  motion 
-22  installation  of  actuators 

-23  if  wanted,  transformation  of  rigid-body  modes  into 
completely  moving  coordinates  of  flight  mechanics 
-24  if  wanted,  further  mode  reduction  for  use  inflight 
simulators 

-25  addition  of  aircraft  pilot  coupling  (APC)  transfer 
functions 

-26  addition  of  electronic  flight  control  system  (EFCS) 

Items  - 1, 2, 3 (partially),  8, 9, 12, 18, 20  are  already  covered  by 
classical  aeroelastic  analyses.  Items  -10, 22, 26  are  to  be  added  at 
least  for  the  analysis  of  ” structural  coupling  of  EFCS”.  The  others 
are  recommended  for  integral  control  design  and  analysis. 

In  the  following  some  of  these  items  are  discussed  in  more  detail. 

2.1  Aerodynamic  Panels  (-3) 

Classical  aeroelastics  most  often  neglects  the  air  loads  of  the 
fuselage.  This  is  not  sufficient  for  an  integral  model.  The  in- 
clusion of  a panel  cross  (s.  Fig.  2.1-1)  with  vertical  and  hori- 
zontal panels  may  be  the  most  simple  approach  if  loads 
caused  by  roll  movements  are  excluded  and  if  all  loads  are 
scaled  (from  panel  to  cylinder  loads).  A scaling  factor  of  0.5 
was  sometimes  found  by  experience. 

2.2  Stiffness  and  Inertia  Corrections  (-6,  -7) 

They  are  second  order  effects,  but  their  realization  is  rather  simple. 
It  simply  requires  the  multiplication  of  the  original  stiffness  and 
inertia  matrices  from  the  right  and  the  left  side  by  a transposition 
matrix  resulting  from  steady  state  structure  deformations. 
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2.3  Control  Modes  (-10) 

A classical  flutter  analysis  includes  about  100  modes  with  the 
lower  frequencies  of  a system  of  about  100000  degrees  of  free- 
dom. Due  to  this  truncation,  local  stiffnesses  may  be  lost.  But  just 
such  local  connection  stiffnesses  may  be  important  if  control  sur- 
faces are  to  be  moved  via  EFCS  by  actuator  loads.  Therefore,  as- 
sumed modes  according  to  Rayleigh-Ritz  are  added.  These  mo- 
des are  defined  by  a unit  deflection  of  the  actuator  and  the 
resulting  deflections  of  the  other  grid  points  of  the  structure.  These 
modes  are  called  control  modes  and  are  to  be  introduced  in  addi- 
tion to  the  classical  rotation  modes  of  the  control  surfaces.  With 
rigid  aircraft  both  modes  would  be  identical. 

More  details  about  the  influence  of  control  modes  are  given  in 
ref.  1. 


2.4  Unsteady  Air  Loads  (—13) 

Until  now  traditional  unsteady  air  loads  are  e.g.  NASTRAN  dou- 
blet lattice  air  loads.  These  are  the  first  harmonics  for  apanel  oscil- 
lating in  an  air  stream.  The  initial  angle  of  attack  is  zero  and  the 
Mach  number  (Ma)  is  constant.  Therefore,  the  following  loads  are 
not  covered: 

- in  plane  loads  (mainly  fore  and  aft  drag  and  side  loads ) 

- unsteady  loads  due  to 

• in  plane  movements  (mainly  fore  and  aft) 

• air  density  Q varying  with  altitude 

• Mach  number  varying  with  altitude  and  velocity. 

For  an  integral  model  this  approximation  is  not  sufficient.  If  better 
air  loads  (e.g.  from  CFD  codes)  are  not  available,  some  correc- 
tions are  necessary.  Details  for  this  are  given  in  ref.  2. 

Here  it  can  be  stated  that  steady  state  air  loads  pO  - also  computed 
by  the  doublet  lattice  method  in  NASTRAN  for  each  aerody- 
namic panel  - can  be  used  to  establish  the  missing  elements.  This 
results  in  the  following: 

For  lift  forces  resulting  from  steady  state  drag  loads 

Apz  = - pOx  * Az/xO  2.4-1 

For  negative  drag  loads  (induced  drag,  neglecting  friction) 

Apx  « (Apz)  • a0  - pOz  • Az/xO  2.4-2 


For  the  other  frequencies  a smoothed  change  to  doublet  lattice  air 
loads  can  be  chosen.  For  a linear  smoothing  this  yields: 

P = Pmeasured  (w)  if  CO  < ©i 
P = Pmeasured  (^)  * ((%—  CO)  / ((Or—  CO i)  + 

+ Pdoublet  lattice  (CO)  * (CO  — tOj)  / CUr  — COi) 
if  COi  < CO  < C0r 
P — Pdoublet  lattice  (w)  if  CO  > (Or 


2.6  Engine  Loads  and  Weight  Forces  (-16,  -17) 

Loads  such  as  gross  thrust,  ram  drag  and  gravity  weight  are  not 
included  in  aeroelastic  analyses,  but  they  are  of  influence  for  flight 
mechanics.  Therefore,  they  must  be  introduced. 

As  to  the  engines,  their  unsteady  load  portion  is  due  to  three  dif- 
ferent reasons: 

First  from  the  unsteady  movement  of  the  engine  position  and  its 
steady  state  loads,  i.e.  due  to  the  change  of  direction  (T’ep0)  and 
due  to  transposition  (V’Eao) 


pO  + Ap  — V'Eo0  ■ T'eEo  ■ p0eE 
resulting  in 

Ap  = - FP(p0eE)  • AroE  2.6-1 

with 


FP(p0eE) 


" 0 0 0 0 -pOz  pOy" 

"Ax 

0 0 0 pOz  o -pOx 

Ay 

0 0 0 — pOy  pOx  o 

Az 

0 -pOz  pOy  o 0 0 

; AroE  = 

Atp 

pOz  o -pOx  0 0 0 

Aft 

-pOy  pOx  o 0 0 0 

Aip 

Second  from  their  derivatives  due  to  velocity  and  altitude  changes 
as  the  air  loads  (2.4—5,  -6) 


^ _ pO(Ma2)  - pO(Mal) 

Ap  ~ Ma2  - Mai 

. / Ax  MaO  . a(z2)  - a(zt) 

y aO  aO  z2  — zi 


Az  + 


+ pQ(q2(z2)  - pQ(e,(zi))  _ Az 
z2  - Zx 
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For  side  loads 

Apy  « -pOz  • A cp  + pOx  • Ay/xQ 
For  in  plane  movements 

p0 

Ap(x)  **  “Ter  * 2 • xO  • Ax 
xUz 
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Third  from  their  derivatives  due  to  control  inputs 


Ap 


p0(n2)-p0(n!) 

n2-nl 
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The  weight  must  be  introduced  if  the  final  equations  are  trans- 
formed from  the  geodetic  (g)  to  the  moving  coordinate  system  (f) 


For  dependencies  from  q and  Ma 


Ap(Q) 


p0(o2(z2))  - p0(Q , (z ! )) 


• Az 


2.4-5 


= gg  * 


-Aft  • cos  cpO  • cos  $0  + Alp  • sin  qpO  ■ cos  00 
Alp  * cos  cpO  * cos  00  - Alp  * sin  00 
-Aqp  • sin  cpO  • cos  00  + A0  * sin  00 


2.6-4 


Ap(Ma)  » 


p0(Ma2)  - pO(Mal) 

Ma2  - Mai 

/ Ax  MaO  a(z2)  - a(zt) 

\ aO  aO  z2  “ zi 
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Another  deficit  to  be  mentioned  here  are  the  doublet  lattice  air 
loads  of  control  surfaces.  Usually,  they  are  too  large  and  should 
be  corrected  to  wind  tunnel  results  of  hinge  moments.  But  the 
corrections  must  also  include  the  airfoil  in  front  of  the  control  sur- 
face. Correction  factors  of  up  to  0.6  are  possible. 


2.5  Smoothing  of  Rigid-Body  Air  Loads  at  Zero  Frequency 
(-15) 

If  steady  state  airloads  from  wind  tunnel  measurements  are  avail- 
able, it  should  be  possible  to  improve  the  doublet  lattice  air  loads 
of  rigid-body  modes  at  zero  frequency.  This  can  be  done  by  sepa- 
rating these  loads  by  an  eigenvalue  analysis.  After  transformation 
of  the  separated  doublet  lattice  air  loads  into  the  coordinate  system 
of  the  measurement,  these  loads  can  be  substituted  by  Fourier  - 
transformed  air  loads  measured  for  zero  frequency  and  the  small- 
est non-zero  frequency. 


2.7  Linearization  and  Transformation  into  Time 
Domain  (-19) 

For  the  purpose  of  controller  design,  the  subsequent  inclusion  of 
nonlinearities  or  the  use  in  flight  simulators,  the  integral  model  is 
required  in  time  domain.  When  starting  with  the  aeroelastics 
model  in  frequency  domain,  a linearization  of  air  loads  and  a 
transformation  are  necessary.  Methods  for  this  were  proposed  by 
Vepa  (s.  ref.  3),  Roger  (s.  ref.  4)  or  Karpel  (s.  ref.  5).  Especially  the 
minimum  state  method  of  Karpel  was  successfully  used,  though 
it  requires  a larger  amount  of  computing.  Fig.  2.7-1  shows  an  ex- 
ample of  good  approximation. 

2.8  Introduction  of  Nonlinear  Rigid-Body  Motion  (-21) 

All  the  improvements  mentioned  above  may  not  be  sufficient  for 
a precise  rigid-body  movement  as  required  by  flight  mechanics. 
Therefore,  a substitution  of  the  improved  ’’aeroelastic”  rigid- 
body  motion  by  the  ’’flight  mechanics”  motion  is  recommended. 
In  this  process  a nonlinear  rigid-body  movement  can  also  be 
introduced. 

Of  course  this  substitution  must  only  be  done  for  the  decoupled 
orthogonal  modes  followed  by  a modal  retransformation  to  the 
original  coordinates. 
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By  starting  from  the  aeroelastics  equation 
xA  = AA  • xA  + BA  • u 
with 

xAr 
xAe 


2.8-1 


xA  = 

r 

tpA  = <j)A-1 


= <J)A 


qAr 

qAe 


_ f XAr  0 
“ I 0 XAe 


1 

qAr 

J 

qAe 

resulting  in  the  orthogonal  coordinates 
qAr 
qAe 

u 

r xpArr  . BAr  + tyAre  • BAe 
+ tpAer  • BAr  + opAee  • BAe 


with  the  eigenvalues 

XAr  (rigid)  and  XAe  (elastic) 

and  defining 

xFr  = (J)Arr  • qAr 

as  intended  by  traditional  flight  mechanics, 

one  obtains  with  equation  2.8-2 
xFr  = (j)Arr  * XAr  • <j)Arr_1  • xFr  + 

+ <})Arr  * (ipArr  • BAr  + tjiAre  • BAe)  • u 

Further  on  it  results  after  some  rearrangements: 


xFr  + 


(J)Are  . XAe  ' ^Aer  c))Are  • XAe  • tpAee 
())Aee  • XAe  * \pAer  <j)Aee  ■ XAe  • t[)Aee 


<J)Are  • (ijiAer  * BAr  + tpAee  ■ BAe) 
4>Aee  • (xJ)Aer  • BAr  + ipAee  • BAe) 


2.8-2 


2.8-3 


2.8-4 


xAr 

"l 

xAe 

4>Aer  • 4)Arr_1 

xAr 

xAe 
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Briefly,  equation  2.8-4  and  -5  can  be  written  as  follows: 


xFr 

xAr 

xAe 


AFrr  0 0 

AFrr  ALrr  ALre 

<t>Ler  * AFrr  ALer  ALee 


xFr 

xAr 

xAe 


BFr 

BFr  + BLr 
(j)Ler  * BFr  + BLe 


xl  = AI  • xl  + BI  ■ u 2.8-6 

That  is  the  integral  equation  of  aeroelastics  and  flight  mechanics. 
If  required,  AFrr  and  BFr  (implying  <j)Frr  and  XFr)  may  be  taken 
from  flight  mechanics  or  the  nonlinear  xFr  = f(xFr,  u)  of  flight 
mechanics  may  be  introduced.  In  any  case,  it  is  the  most  reason- 
able compromise.  It  keeps  the  eigenvalues  of  the  rigid  a/c  Xr  of  the 
discipline  defining  xFr  and  keeps  the  eigenvalues  of  the  elastic 
structure  Xe  of  aeroelastics  unchanged  as  can  easily  be  proven. 
Even  simulations  in  a moving  simulator  with  elasticity  included 
can  be  performed.  The  eigenvector  matrix  of  AI  reads: 

^(4>Frr  or  <|)Arr)  0 0 

(cj)Frr  or  tpArr)  -1  (|)Are 

<|)Aee  • c|)Are_1  * (<j)Frr  or  oJiArr)  -<t>Aer  • <t>Arr  cJ>Aee 

and  their  eigenvalue  matrix: 


(XFr  or  XAr) 

0 

XAe 
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2.9  Installation  of  Actuators  (-22) 

Actuators  are  highly  nonlinear  elements,  mainly  due  to  the  vis- 
cous damping  effect  of  the  oil  flow,  the  overlapping  of  the  throttle 


orifices  and  the  quantization  effects  of  the  digital  servo  control 
loop.  Therefore,  actuators  have  to  be  substituted  by  their  fre- 
quency response  functions  (in  the  following  called  transfer  func- 
tions), if  linear  systems  are  to  be  used.  This  includes  the  first  har- 
monic but  neglects  all  higher  orders  or  subharmonic  responses. 
Then  a decision  has  to  be  taken  on  the  actuator  model.  There  are 
two  possibilities  for  actuators  with  closed  control  loop. 

The  stroke  model: 

zA  = H2A,u  • u 2.9-1 

and  the  load  model: 

zA  = HzAu  * u + HzApA  • pA  2.9-2 

where  zA  = realized  actuator  stroke 

u = commanded  actuator  stroke 
pA  = actuator  action  load 
H. , . = transfer  function 

(Hza,u  may  be  called  ’’frequency  response  function 
of  control”  and  Hza,pa  ’’frequency  response  function 
of  disturbance”) 

and  after  some  rearrangements  of  equation  2.9-2  for  the 
loads  model: 

pA  = pH  - kA  • zA  2.9-3 

with 


pH  — HpHu  ■ u + HpH  zA  • zA 

HpH.u  "HzAjpA 

Hph,za  = HzApA  + kA 

kA  = chosen  more  or  less  arbitrarily 

Both  models  have  their  advantages  and  disadvantages. 

2.9- 1  The  stroke  model 

This  is  the  simpler  model  and  the  one  most  often  used.  It  requires 
the  elimination  of  za>  the  actuator  stroke  or  deflection,  from  air- 
craft structure  equations  and  its  transfer  to  the  right  side  of  these 
equations  to  make  it  available  as  excitation. 

The  remaining  homogenous  equations  of  structure  on  the  left  side 
do  no  longer  include  actuators.  They  are  only  valid  for  rigid  actua- 
tors or  rigid  connections  between  control  surface  and  aircraft. 
Therefore,  two  eigenvibration  calculations  are  required.  The 
usual  one  for  classical  flutter  analyses  without  electronic  control 
but  with  flexible  actuators  and  the  other  one  for  flutter  analyses 
with  electronic  control  and  with  rigid  actuators.  Most  probably, 
the  eigenfrequencies  and  mode  shapes  of  these  two  models  will 
vary  to  a greater  extend. 

The  models  represent  two  different  working  points.  The  one  with 
rigid  actuator  may  be  too  far  away  from  the  conditions  needed 
later  when  the  actuator  is  flexible  and  moved  by  electronic  con- 
trol. 

Besides,  it  should  also  be  mentioned  that  in  case  of  polynomial 
approximations  of  the  transfer  function  of  the  stroke  model,  the 
order  of  the  denominator  polynom  must  be  two  orders  higher  than 
the  nominator  polynom  to  get  accelerations  of  actuator  stroke 
from  commanded  stroke. 

2.9- 2  The  loads  model 

This  model  is  the  more  complicated  one.  Its  advantage  is  that  it  al- 
lows the  inclusion  of  actuator  flexibility  existing  under  dynamic 
conditions.  The  model  operates  as  areal  linear  spring  andadistur- 
bance  load  pH  in  parallel.  If  the  spring  term  -kA  * za  of  equation 

2.9- 3  is  brought  to  the  left  side  of  aircraft  structure  equations,  a 
classical  eigenvibration  calculation  with  flexible  actuators  can  be 
performed  and  one  set  of  modes  can  be  used  for  flutter  analyses 
with  and  without  electronic  control.  This  is  a second  advantage  of 
the  loads  model.  The  stiffness  kA  ( and  perhaps  even  an  additional 
damping  term)  can  be  taken  from  ground  vibration  tests.  A third 
advantage  could  be  the  fact  that  the  introduction  of  Hph,z a in  prin- 
ciple allows  a stiffness  correction  even  after  modal  truncation,  as 
explained  already  in  ref.  1. 

The  control  term  Hph,u  * u and  the  difference  between  complex 
and  real  negative  actuator  spring  (HzApA  + kA)  • zA  remain  on 
the  right  side  of  aircraft  equations.  An  example  of  the  negative 
complex  spring ’’stiffness”  HnA,zAts  given  in  Fig.  2.9-1,  together 
with  a chosen  real  substitute  4ca.  The  differences  are  not  negligi- 
ble. 

Fig.  2.9-2  gives  some  information  about  the  influence  of  different 
terms  of  the  actuator  transfer  function.  It  shows  a vertical  wing  tip 
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acceleration  response  due  to  elevator  movement  in  flight.  Three 
cases  are  shown.  The  first  one  for  the  complete  actuator  model 
with 

^pH,u  — ^zA,pA  ^zA,u 
^pH,zA  ^zA,pA  ^A 

the  second  one  arbitrarily  with 
HpH,u  = ~HzApA  ■ HzAu 

Hph,zA  - 0 

the  third  simplified  one  with 
HpH,u  = ‘ H2a,u 

HpH,zA  = 0 

In  the  figure  the  peaks  of  structural  response  are  clearly  visible. 
The  second  case  shows  large  deviations  at  about  17  Hz,  while  for 
the  third  case,  only  small  deviations  at  about  8 Hz  are  visible. 
So  it  can  be  concluded  that  the  complex  frequency  response  func- 
tion of  disturbance  Hza,pa  has  a larger  influence  and  should  be 
precise  enough,  if  used. 

2.9-3  Scatter 

Another  important  factor  of  actuator  models  is  the  large  number 
of  influence  parameters,  such  as 
Q temperature,  environment 
Q backlash,  wearing 
Q supply  and  return  pressure 
Q servo  control  law  (gain,  sample  rate) 

□ activated  mode 

□ amount  of  mass  to  be  moved 

□ static  load  in  working  point 

Q static  position  in  working  point 

□ amplitude  and  frequency  of  command 

□ amplitude  and  frequency  of  external  load 

□ failure  conditions 

• supply  pressure  decrease 

• leakage 

• wrong  activation  status 

• desynchronization  of  multiple  actuators 

• mechanical  blocking  or  fracture 

• air  in  hydraulic  oil 

Therefore,  a broad  scatter  band  of  transfer  functions  is  to  be  ex- 
pected. The  influence  of  this  scatter  is  shown  in  Fig.  2.9-3.  The 
same  transfer  function  as  the  one  shown  in  Fig.  2.9-2  is  given  for 
the  complete  actuator  model,  together  with  its  upper  and  lower 
tolerance  (mainly  due  to  oil  temperatures  of  90°,  35°  and -15°). 
The  differences  are  dramatic.  As  a consequence,  stability  and  re- 
sponse analyses  must  include  a reasonable  amount  of  scatter. 

2.10  Transformation  into  Completely  Moving 
Coordinates  (-23) 

Flight  mechanics  use  a coordinate  system  which  moves  com- 
pletely with  deflection.  Aeroelastics  use  an  inertial  geodetic  coor- 
dinate system.  The  integral  model  derived  from  aeroelastics  starts 
with  the  inertial  system.  Therefore,  a transformation  may  some- 
times be  necessary.  For  the  sake  of  completeness,  it  should  be 
mentioned  here  that  any  transformation  into  the  flight  mechanics 
system  adds  so-called  Euler  terms  to  the  equations.  They  can  sim- 
ply be  derived  as  follows: 


with 

sgt  velocity  with  components  measured  in  the 
inertial  coordinate  system  g 
sft  same  velocity,  but  its  components  are  now 
measured  in  the  completely  moving  coordinate 
system  f. 

TfAg  matrix  of  coordinate  transformation  from  g to 
f with  Euler  angles 

the  following  results: 

^ft  — ^fAg  ^gt  + TfAg  • Sgt  2.10—2 

with  the  first  term  being  the  Euler  term  (derivated  without  cross 
products  or  special  derivation  symbols  for  derivation  in  moving 
coordinate  systems!).  The  Euler  term  is  not  negligible  since  sgtin- 


cludes  the  velocities  of  the  aircraft.  All  the  other  terms  contain 
small  variables. 

2.11  Addition  of  Aircraft  Pilot  Coupling  (APC),  (-25) 

In  classical  aeroelastics  the  movement  of  an  aircraft  is  determined 
exclusively  by  the  aircraft  structure  and  the  air.  But  with  highly 
flexible  a/c,  a pilot  in  the  loop  must  also  be  investigated. 

If  e.g.  the  floor  of  the  aircraft  cabin  is  vibrating,  this  vibration  is 
transferred  via  seat,  pilot  body  and  arm  to  the  side  stick.  The  side 
stick  is  unintentionally  moved  and  introduces  signals  to  the  actua- 
tors deflecting  the  control  surfaces  which  excite  the  aircraft  via  air 
loads  and  closes  the  loop  of  vibration  (s.  Fig.  2.11-1). 

This  is  not  only  of  theoretical  interest,  critical  cases  - called  PIO 
cases  before  - are  known  from  several  aircraft. 

Therefore,  the  complete  integral  model  should  include  this  bio- 
mechanical closed  loop  circuit  of  APC. 

This  can  be  done  if  a transfer  function  from  aircraft  cabin  floor 
vibration  to  stick  input  is  available.  Such  functions  were  recently 
measured  and  some  are  published  in  ref.  6.  Fig.  2.1 1-2  shows  an 
example. 

If  everything  mentioned  up  to  now  is  introduced,  the  integral 
model  is  ready  for  design  and  application  of  integral  control. 

2.12  Results  from  an  Integral  Model 

An  aircraft  model  including  most  of  the  effects  mentioned  above 
was  established  and  its  characteristics  were  studied. 

Based  on  flight  mechanics,  three  different  sets  of  linearized  rigid- 
body  models  were  available.  An  early  one  FI Lr,  an  improved  one 
with  ’’rigid”  aerodynamic  coefficients  F2Lr  and  one  with  ’’elasti- 
fied”  aerodynamic  coefficients  F2Le. 

Based  on  aeroelastics,  three  different  sets  of  models  are  available. 
The  first  one  A(qr)  includes  only  rigid-body  modes,  the  second 
one  A(qs)  includes  the  first  six  symmetric  and  six  antisymmetric 
elastic  modes,  and  the  third  one  A(qc)  includes  two  times  27 
instead  of  12  elastic  modes. 

2. 12- 1  Eigenvalues 

Table  2. 12-1  gives  a survey  of  the  eigenvalues  (frequency  fn  and 
damping  £n)  of  the  most  important  rigid-body  modes. 

For  the  ’’short  period”  the  frequencies  of  flight  mechanics  are  in 
general  smaller  than  the  frequencies  of  aeroelastics,  while  the 
dampings  show  opposite  tendencies.  Elastic  a/c  show  smallerfre- 
quencies  and  larger  dampings  for  both  disciplines.  The  influence 
of  elasticity  is  not  negligible. 

For  the  ’’dutch  roll”  the  frequencies  are  in  better  agreement,  where 
as  the  dampings  of  flight  mechanics  show  larger  scatter.  The  influ- 
ence of  elasticity  is  smaller. 

The  ’’phygoid”  is  rather  sensitive  and  a larger  scatter  is  visible  for 
the  values  of  flight  mechanics.  Nevertheless,  the  frequencies  of 
aeroelastics  do  not  differ  so  much.  The  influence  of  elasticity  is 
visible. 

The  other  aperiodic  rigid-body  modes  show  relatively  good 
approximation.  The  increase  of  the  number  of  elastic  modes  to  54 
has  minor  influences  on  the  frequencies  and  dampings  of  the  first 
12  elastic  modes. 

2.12— 2  Response  to  Control  Inputs 

Figures  2.12-1  to  -4  show  examples  of  time  response  of  rigid- 
body  states 

• pitch  acceleration  due  to  elevator  ’’step”  (+1  degree) 

• roll  acceleration  due  to  aileron  ’’stair”  (+l°,0°,-lo,0o) 

• yaw  acceleration  due  to  rudder  step  (+1°) 

Figure  2. 12-1  compares  the  different  models  of  flight  mechanics 
with  the  reference  model  of  aeroelastics  A(qs).  The  differences, 
mainly  in  yaw  damping,  and  the  influence  of  elasticity  are  clearly 
visible. 

Fig.  2.12-2  compares  the  influence  of  a different  number  of  elas- 
tic modes.  It  shows  the  influence  of  elasticity  more  clearly.  Espe- 
cially with  pitch  and  roll  acceleration,  it  is  visible  that  with  elastic 
modes  the  vibrations  move  around  their  own  average  curve  but 
not  around  the  curve  of  the  simple  rigid  model.  This  proves  again 
the  influence  of  elasticity  on  rigid-body  movement. 

Fig.  2.12-3  compares  the  results  of  the  integral  model  I(A(qs)), 
i.e.  (s.  § 2.8)  rigid-body  behavior  determined  by  aeroelastics,  xFr 
= xAr(A(qs))  with  I(F2Lr),  i.e.  rigid-body  behavior  determined 
by  flight  mechanics,  xFr  = xF(F2Lr). 

The  differences  are  rather  small.  They  result  from  the  different 
rigid-body  models. 

In  addition,  one  can  see  clearly  that  the  same  elastic  vibration  is 
superimposed  to  the  applied  rigid-body  movement.  This  proves 
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again  that  elasticity  can  be  superimposed  by  applying  equation 
2.8-5  if  a good  rigid-body  model  is  available. 

Finally  Fig.  2.12-4  shows  the  measurable  acceleration  at  the  pilot 
seat  (model  (A(qs))  due  to  control  commands.  Here  the  influence 
of  elasticity  is  rather  large  as  already  shown.  The  response  ampli- 
tudes due  to  rudder  step  are  nearly  twice  as  large  as  those  of  rigid 
a/c,  the  same  applies  to  the  elevator  step  which  intro  duces  a strong 
response  peak  right  at  the  beginning. 

Naturally,  such  a large  elastic  a/c  response  would  be  reduced  in 
reality  by  electronic  filtering  or  by  mechanical  reduction  due  to 
slow  actuators.  In  this  simulation  no  reduction  was  introduced. 
The  transfer  functions  of  actuators  were  set  to  1 =constant  with  the 
intention  to  show  the  potential  of  the  aircraft. 

3.  THE  INTEGRAL  CONTROLLER 

3.1  Objectives 

There  are  three  different  tasks  of  the  integral  controller  to  be  dis- 
tinguished: flight  control,  load  control  and  mode  control.  Their 
objectives  are  rather  different. 

In  general  one  can  name: 
flight  control: 

- flight  stability 

- aircraft  protection 

- maneuverability 

- defined  flight  path  following 

- reduction  of  workload  of  pilot 
load  control: 

- reduction  of  limit  loads 

- reduction  of  fatigue  loads 

- avoidance  of  critical  failure  conditions 
mode  (or  structural  vibrations)  control 

- aeroelastic  stability 

- minimization  of  structural  vibrations 

- improvement  of  passenger  comfort 

All  these  objectives  interfere  with  each  other,  some  are  even  con- 
tradictory to  each  other,  and  most  of  them  can  be  subdivided  into 
a larger  number  of  subobjectives. 

Nevertheless,  they  can  be  detailed  and  - with  some  exceptions  - 
given  as  criteria  in  a mathematical  formulation.  So  an  integral  con- 
troller can  be  designed  in  a multicriteria  optimization  process, 
with  targets  and  constraints. 

For  the  purpose  of  mode  control,  which  is  of  main  interest  in  this 
context,  one  can  distinguish  three  different  types  of  controller  ob- 
jectives: flutter  suppression,  flutter  margin  augmentation  and  re- 
duction of  vibration  level  for  comfort.  Flutter  suppression  would 
cure  a flutter  case  within  the  flight  envelope.  Flutter  margin  aug- 
mentation would  cure  aflutter  case  in  an  area  for  which  the  aircraft 
is  not  designed  but  which  must  be  covered  as  required  by  law. 
Comfort  control  would  reduce  the  natural  structural  vibrations 
caused  by  air  turbulence  down  to  a comfortable  level  on  a natu- 
rally stable  aircraft. 

All  mode  controllers  aim  at  an  increase  of  the  damping  of  structure 
modes.  Their  main  difference  is  the  required  level  of  safety. 

3.2  Requirements  and  Restrictions  for  Mode  Control 

3.2- 7  Airworthiness  Requirements 

For  mode  controllers  the  main  safety  requirements  are  also  air- 
worthiness requirements.  For  civil  aircraft  and  structural  aspects, 
these  are  the  FAR  or  JAR  25  requirements,  here  § 25.629  ’’Flutter, 
Deformation  and  Fail  Safe  Criteria”  and  § 25.302  ’’Interaction  at 
Systems  and  Structure”.  If  one  looks  at  the  flutter  speed  to  be 
reached  a literal  interpretation  of  the  requirements  leads  to  Fig. 

3.2- 1  for  its  definition.  The  required  speed  depends  on  the  sys- 
tem’s probability  of  being  in  failure.  So  one  can  read  in  Fig.  3.2-1 
that: 

Q flutter  suppression  controllers  need  a probability  of  being  in 
failure  state  of  <10-9,  since  they  have  to  shift  the  flutter 
speed  from  below  VD  up  to  1.15  VD 
Q flutter  margin  augmentation  controllers  need  < 1 0-5 , since 
they  have  to  shift  from  above  VD  up  to  1.15  VD 
Q comfort  controllers  need  <1,  since  here  the  natural  a/c 
reaches  already  1.15  VD.  If  the  system  fails,  there  is  no  safety 
problem  but  only  a comfort  problem  which  is  not  the 
subject  of  these  paragraphs. 

A special  situation  arises  if  a high  speed  protection  (HSP)  system 
and  a mode  control  (MC)  system  are  used  in  combination  to  keep 
V < Vc  and  to  avoid  flutter  within  Vc<Vp<Vo.  If  both  systems 
have  an  independent  probability  of  being  in  failure  state  of  <10-5, 


in  combination  they  can  reach  a value  of  < 10~9,  as  required.  But 
nevertheless,  Vp<Vbif  the  MC  system  fails,  although  VF  cannot 
be  reached  due  to  the  HSP  system. 

3.2-2  Mode  Control  Design  Criteria 

Any  controller  should  be  safe,  efficient,  robust  and  should  not  dis- 
turb other  positive  features  of  the  controlled  system. 

For  safety  reasons,  the  a.m.  JAR  requirements  definedby  the  flut- 
ter speed  to  be  reached  must  be  fulfilled.  But  there  are  other  safety 
criteria  too. 

First  the  amount  of  modal  damping  which  indicates  a certain  sta- 
bility robustness  against  scatter,  ageing  or  disturbances. 
Second  the  phase  and  gain  margins  (well  known  for  single  input/ 
single  output  (SISO)  - systems)  which  show  the  robustness 
against  deviations  from  nominal  conditions  or  against  minor  fail- 
ures not  especially  mentioned. 

The  efficiency  or  performance  of  the  controller  can  very  simply 
be  defined  inmost  cases  by  the  flutter  velocity  or  the  modal  damp- 
ing to  be  reached. 

The  robustness  is  a more  complex  issue.  One  aspect,  the  ’’safety 
robustness”  was  already  mentioned  before.  Another  aspect  is 
’’performance  robustness”  against  disturbances  or  against  opera- 
tion conditions.  The  number  of  the  latter  one  is  huge.  It  includes, 
among  others,  a/c  velocity,  altitude,  flap/sl at  setting,  payload,  fuel 
and  their  distributions.  Each  combination  of  these  conditions  de- 
fines a load  case  to  be  included  in  the  controller  design. 

Finally  it  must  be  shown  by  the  evaluation  of  time  response  due 
to  step  input  or  frequency  response  due  to  random  excitation  that 
the  controller  neither  introduces  unwanted  features  nor  reduces 
wanted  features,  such  as  good  a/c  maneuverability. 

In  the  following  some  requirements  are  detailed  as  a first  proposal 
for  a multi-input  / multi-output  (MIMO)  controller. 

□ Flutter  control: 

For  modal  damping  ”§n”  of  each  of  n eigenvalues  shall  be 
^n  > 0.5%  (%  of  critical)  and 

^ > -0.01  [%  / kts  CAS] 

This  guarantees  stability. 

□ Comfort  control: 

The  maximal  singular  value  ”0”  of  a number  m of  pick-up 
signals  ”ym”  due  to  gusts  ”w”  in  frequency  range  ”fi” 
to  ”f2”  shall  be  minimal  and  at  least  smaller  than  a limit 
value  Guym.w  ie. 

max  (a  (Hym,w  (ioo)))  ->  min 
fl-f2 

max  (0  (Hym,w  (ito)))  < GHym,w 

0-s-  00 

where  Hym,w  (ito)  is  the  transfer  function 
from  gust  w to  the  pick-up  signal  ym. 

This  should  reduce  the  response  to  turbulence. 

□ Robustness  of  control  (based  on  ref.  7): 

The  phase  margin  of  the  open  loop  transfer  function 
Hu,u  (ito)  of  all  actuator  signals  u shall  be 

^ T80  K'  “ I ^ ’ arc™  (n™  fe  (1  + Hu,u  (ito)))  /2)  | 

the  gain  margin  of  the  open  loop  transfer  function  Hu,u  (s) 
of  all  actuator  signals  u shall  be 
1 / (1-min  (0  (1  4-  Hu,u  (ito))))  > 2 
and 

1/(14-  min  (0(14-  Hu,  u (ito))))  < 0.66 
where  0 is  the  minimum  singular  value. 

This  should  ensure  safety  robustness. 

If  these  criteria  were  really  met  in  the  whole  envelope  of  velocity 
and  altitude  as  well  as  in  all  loading  conditions  of  fuel  and  pay- 
load,  a good  controller  should  be  reached. 

But  these  requirements  are  probably  too  stringent  and  some  moni- 
tored exceptions  are  therefore  necessary.  At  the  moment,  there  is 
insufficient  experience  available  to  present  more  details  on  this  is- 
sue. But  one  thing  is  beyond  question:  some  requirements  of  this 
type  are  necessary. 

The  mere  definition  of  the  required  flutter  speed,  as  given  in  FAR 
25.302,  is  not  sufficient.  At  least  a certain  amount  of  robustness 
is  necessary  for  safety  to  cover  smaller  scatter  or  deviations  from 
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nominal  conditions  which  are  not  defined  as  separate  failure 
states. 

3. 2- 3  Reliability  / Failure  Probability 

A good  reliability  or  low  failure  probability  can  be  reached  by  dif- 
ferent methods.  For  example,  redundancy  of  hardware  installa- 
tions and  signal  monitoring  are  well  known.  For  MIMO  control- 
lers studied  here,  special  attention  should  also  be  paid  to  the 
minimization  of  interactions  between  the  different  signal  lanes  of 
the  controller. 

3.2- 4  Safety  Restrictions  for  Mode  Control 
This  could  include: 

□ Low  pass  filter  for  the  flight  control  signal  part  and  band  pass 
filter  for  the  signal  part  of  mode  control.  This  prevents  a 
failure  in  the  more  complicated  flight  controller  part  from 
causing  flutter  problems  and  a failure  in  the  mode  controller 
part  from  influencing  the  flight  path. 

Q Decoupling  network  with  diagonal  structure  of  the  controller 
transfer  matrix  of  MIMO  systems.  This  could  e.g.  be  done  by 
superposition  of  different  effective  SISO  systems  which  do 
not  interfere  with  each  other. 

If  all  SISO  systems  are  stable  and  effective  for  themselves,  a 
failure  and  a switch-off  of  one  of  these  lanes  will  not  lead  to 
a complete  loss  of  the  controller. 

Admittedly,  such  measures  cannot  always  be  taken.  Especially  a 
separation  by  filter  does  not  work  if  rigid-body  frequencies  and 
elastic  mode  frequencies  are  close  together.  But  if  it  works,  these 
filters  should  be  used  under  all  circumstances.  Missing  filters 
present  an  unnecessary  risk. 

3.3  Design 

The  design  of  the  integral  controller  is  a multidisciplinary  process 
of  multicriteria  optimization.  Fig.  3.3-1  may  give  a survey  of  its 
principle  structure. 

First  an  integral  model  of  the  a/c  must  be  established  and  made 
available  in  linear  and  nonlinear  form  for  all  loading  and  flight 
conditions  to  be  covered. 

Second  a complete  list  of  criteria  is  to  be  established,  including  all 
objectives  of  the  different  disciplines. 

Then  the  question  arises  how  the  controller  should  be  structured. 
In  most  cases  experience  answers  this  question,  otherwise, 
theoretical  methods  of  robust  optimal  control  law  design  as  e.g. 
the  Hqo  method  are  available  to  find  a good  structure. 

If  the  structure  is  known  its  parameters  can  be  optimized  by  a 
multiobjective  optimization  strategy. 

Finally  the  resulting  controller  must  undergo  a special  assessment 
by  the  different  disciplines,  including  simulations  with  pilot  in  the 
loop. 

Of  course,  one  has  to  keep  in  mind  that  by  increasing  the  number 
of  aims  the  efficiency  of  the  controller  will  be  reduced. 

It  is  still  impossible  to  get  everything  simultaneously! 
Modifications  of  controller  structure,  of  criteria  or  even  of  control 
surfaces  or  systems  may  be  necessary. 

The  first  three  iteration  steps  of  an  example  of  such  a multicriteria 
optimization  is  shown  in  Fig.  3.3-2  (s.  ref.  8)  where  a limited 
number  of  criteria  was  studied.  One  can  see  that  finally  all  criteria 
are  met.  The  program  stops  when  ’’pareto  optimal”  results  are 
reached,  i.e.  when  a value  of  one  criteria  can  only  be  improved  by 
impairing  another  one. 


4.  CONCLUSION 

□ Increasing  a/c  flexibility  requires  integral  controllers  for 
flight,  load  and  mode  control. 

□ The  controller  design  must  be  based  on  an  integral  a/c  model 
which  covers  the  requirements  of  flight  mechanics,  loads  and 
aeroelastics. 

Q The  integral  model  can  be  based  on  a classical  aeroelastics 
model  if  several  improvements  are  introduced. 

□ Linear  or  nonlinear  rigid-body  motion  of  traditional  flight 
mechanics  can  be  coupled  with  structural  vibrations  of 
traditional  aeroelastics  without  restrictions. 

□ APC  can  be  included. 

Q Controller  design  should  be  performed  by  a multiobjective 
parameter  optimization  to  reach  an  optimum  for  all  criteria. 

□ Safety  and  design  criteria  for  mode  control  need  further 
discussion. 
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FIG.  2.1-1:  MODELING  OF  FUSELAGE  FOR  AIR  LOADS 
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FIG.  2.7-1  TRANSFER  FUNCTION  FROM  ELEVATOR 
TO  FRONT  FUSELAGE  ACCELERATION 


FIG.  2.9-1:  ACTUATOR  TRANSFER  FUNCTIONS 
OF  DISTURBANCE  OR  FLEXIBILITY 
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FIG.  2.9.— 2 
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FIG.  2.9-3 
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INFLUENCE  OF  ACTUATOR  SCATTER 
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FIG.  2.11-1:  AIRCRAFT  PILOT  COUPLING  MECHANISM 
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FIG.  2.1 1 -2:  TRANSFER  FUNCTION  FROM  A/C  FLOOR  VIBRATION 

TO  SIDE  STICK  MOVEMENT 


[••/Oop]  pp  m [••/6«P]  PP  md  [••/flap]  PP  |«d 


9-9 


phygoid 


0.01076 

5.22 


0.01160 

7.30 


0.01232 

8.14 


dutch 

roll 


0.1596 

10.70 


0.1551 

10.87 


0.1480 

10.01 


short 

period 


0.2442 

48.76 


0.2202 

48.96 


0.1876 

51.87 


0.007706 

0.1640 

0.1694 

-8.05 

11.20 

55.08 

0.009915 

0.1685 

0.2105 

2.65 

6.83  ' 

i 

52.72 

0.01140 

0.1630 

0.1391 

2.01 

2.17 

68.26 

TABLE  2.12-1:  EIGENVALUES  OF  INTEGRAL  MODELS 
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FIG.  2.12-1:  ACCELERATION  RESPONSE  TO  CONTROL  COMMAND 

(Comparison  of  Data) 
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FIG.  2.12-4:  PILOT  ACCELERATION  IN  ELASTIC  A/C 


1.15  • (VDp;  (or  1.15-  VCp;  or  V(MD);  or  V(Mc+0.5))) 


Probability  of  being  in  failure  state  Qj  = pj  *Tj  with 
average  time  Tj  spent  in  failure  condition  and  pj 
probability  of  occurrance  of  failure  mode  per  flight  hour. 

e.g.: 

VDu  - VCu  « 60  kCAS  unprotected 
VDp  - VCp  « 35  kCAS  protected 


FIG:  3.2-1:  REQUIRED  FLUTTER  SPEED  (NPA  25C-199/ACI  25.302  § 4.1. 2.2) 
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